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We reexamine large CP-violating phases in the general Minimal Supersymmetric Standard Model, 
as well as more restricted models. We perform a detailed scan over parameter space to find 



Oh' 

solutions which satisfy the current experimental limits on the electric dipole moments of the 
electron, neutron and 199 Hg atom, exploring the allowed configurations of phases and masses, and 
we attempt to quantify the level of tuning of the parameters necessary to populate the regions 
of cancellations. We then consider the measurement of CP-violating phases at a future linear 
collider. We find that measurements of chargino and neutralino masses and production cross- 
sections allow for a determination of 0i(the phase of Mi) to a precision of 7r/30, while the EDM 
constraints require that 9^ be too small to be measured. Using the EDM constraints we find that 
the CP-even model parameters and the phase <pi can be determined at a Linear Collider with 
400 GeV cm. energy. As long as some information on the size of \fi\ is included in the observables, 
a measurement of <pi is guaranteed for 0i > 7r/10. To unambiguously identify CP violation, we 
construct CP-odd kinematical variables at a linear collider. However, the CP asymmetries are 
rather small, typically about 0.1 — 1.5%, and it will be challenging to experimentally observe the 
predicted asymmetries. 
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1 Introduction 



In the last two years, there has been a resurgence of interest in CP-violating phases in supersym- 
metric theories. Although in many cases the size of the CP- violating phases is strongly constrained 
by limits on the electric dipole moments (EDMs) of the electron, neutron and 199 Hg atom, possi- 
ble cancellations between different contributions to the EDMs can significantly weaken the upper 
bounds on the phases |J, even potentially allowing phases of 0(1) P, §, typically in small, often 
finely tuned regions. This observation has led to an explosion of papers exploring the consequences 
of the presence of large phases, many papers finding some regions of parameter space where large 
phases are of consequence to the system under study, but few papers simultaneously imposing the 
(still severe) constraints from EDMs. One of the most interesting questions is whether one will 
be able to measure the phases in chargino and neutralino production at a future linear collider 
0, Even if one expects that weakly interacting supersymmetric partners are going to be found 
at the Tevatron and at the LHC || [F], the model parameters have to be determined in detail in a 
precision experiment. A future high luminosity Linear Collider is the optimal tool to extend the 
LEP Standard Model precision measurements to other models. The CP-conserving set of under- 
lying MSSM parameters can be easily determined from physical masses alone, cross sections and 
asymmetries H . The effects of CP- violating phases on neutralino and chargino observables can be 
considerably smaller and therefore require a more ambitious collider energy and luminosity [[|, |JJ . 

In this paper we reexamine this question, taking into account the limits from EDM experiments, 
including the recent improved limits on the 199 Hg EDM [[|. We perform a detailed scan of 
parameter space in the general (non-unified) MSSM, including phases and masses, to find the 
regions satisfying the EDM constraints. Using these constraints we show how the phase <pi can 
be extracted from a set of CP-even variables (masses and total cross sections) at a future Linear 
Collider. Compared to the previous analysis Q we can reduce the requirements on the energy 
as well as on the luminosity of the Linear Collider. Moreover, we show how one could use an 
extended set of parameters, which allows for an independent extraction of and further relax 
the requirement on the Collider design parameters. Finally, we construct CP-odd variables at 
e + e~ linear colliders. We find that the CP asymmetry is typically about 0.1 — 1.5%. 



2 CP Violation and Electric Dipole Moments 

The Minimal Supersymmetric Standard Model (MSSM) contains several sources of CP violation 
not present in the Standard Model. In the most general flavor non-preserving SUSY model, there 
are over 40 new complex phases flO| , although most of the new parameters are very strongly 
constrained by limits on flavor violating processes. We assume a Peccei-Quinn symmetry |T] 
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which sets the coefficient 9 of the GG term to zero, up to small corrections coming from higher 
dimensional CP- violating operators, which shift the minimum of the axion potential |12|]. In the 



MSSM, new complex phases arise in the Higgs mixing mass /i, as well as in the soft SUSY-breaking 
terms in the Lagrangian: the trilinear scalar mixing masses A iy the bilinear Higgs mixing parameter 
B, and the three gaugino mass parameters Mj. Not all of these phases are physical, however, and 
some or most may be removed by field redefinitions, depending on the model. In mSUGRA, for 
example, the phases 0, of all the Mi may be removed by an R-rotationQ at the unification scale 
Mx-i the Ai are united to a common A Q at Mx, and a rotation of one of the Higgs fields can be 
used to set 6^ + Ob = 0, so that the Higgs vevs are real. This leaves only two physical phases in 
mSUGRA, which can be taken to be 6U and 6^. 

In more general models, where the Mj do not unify, only one of the gaugino phases can be 
rotated away, which we take to be the phase of M 2 . In our more general analyses, we will also 
take independent phases in the trilinear parameters A e , A u , Ad and A t . This leaves us with seven 
physical phases in our most general model. We emphasize that these are the phases present after 
phases have been removed in the field redefinitions described above. Therefore when we constrain 
below the phase of the \x parameter to be very small, it is in fact some combination of the phases 
in the original parameterization which is restricted. This is particularly important in models 
where there is a correlation between the phase of B and the other phases in the model, such as in 
mSUGRA, as we discuss further below. 

It is well known that the additional sources of CP violation in supersymmetric models can 
contribute to the electric dipole moments of the neutron and electron JT3], II and mercury atom 
TEfl . The very tight experimental limits on these quantities |17], || 



d n < 1.1 x KT 25 e-cm (1) 
4 < 4 x KT 27 e • cm (2) 
d Ug < 2 x KT 28 e • cm (3) 

then impose severe constraints on the CP-violating phases in SUSY models. To suppress the 

EDMs, either large scalar masses ( > 1 TeV) or small phases (of the order 10 -3 , when all SUSY 

masses are of order 100 GeV) are typically required. However, as pointed out in |T], |T%fl , such 

large scalar masses are cosmologically problematic, and the addition of cosmological constraints 

to the mix implies that (some of) the CP-violating phases are constrained to be quite small when 

the LSP is a dominantly S-type neutralino. 

The electron EDM receives contributions from chargino and neutralino exchange diagrams, 

shown in Fig. [l]. The full expressions for the dependence of the induced electron EDM on the 

1 Under an R-rotation of the fields, the Ai, B and Mi effectively pick up a common phase, while the other soft 
SUSY breaking masses and jj, are left unchanged. 
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SUSY masses and phases can be found in [|14], |2|]. The dominant contribution is typically from 
the chargino exchange diagram and is proportional to sin 6^. Thus the primary constraint coming 
from the electron EDM limits is an upper limit on the phase of fi. The subdominant neutralino 




Figure 1: Diagrams contributing to the electric dipole moment of the electron. The SUSY phases 
enter at the crossed vertices. 



exchange piece has a more complicated dependence on the SUSY phases, including both pieces 
proportional to sin^ and shi7 e , where 7 e = arg(A e + /i* tan/3). Cancellations can occur between 
the neutralino chargino exchange contributions, and this serves to weaken the absolute limits on 
the SUSY phases, although 9 ^ of 0(1) requires either severe fine-tuning of parameters or a very 
heavy spectrum, as we will see below. 

The neutron EDM is considerably more complicated, and until recently, the computation 



14| . [2j of the neutron EDM induced by SUSY phases has been plagued by very large theoretical 
uncertainties. The SUSY phases contribute both to the EDMs and color EDMs (cEDMs) of the 
quarks, and in the last year, the contribution to the neutron EDM both from the induced 9 due 



to the color EDMs of the quarks |20| and from the quark EDMs and cEDMs themselves |TL| have 
been reliably calculated using QCD sum rules, allowing a reduction in the theoretical uncertainty. 
The overall effect is to reduce slightly the predicted neutron EDM, and with smaller error bars. 
The quark EDMs receive contributions from the chargino and neutralino exchange diagrams of 

\ 

Figure 2: Diagram contributing to the electric dipole moment of the down quark. The SUSY 
phases enter at the crossed vertices. 



Fig. [I] (where the photon can now also connect to the sfermion line in the chargino diagram), 
as well as from the gluino exchange diagram of Fig. 0. The gluino exchange contribution to the 
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d-quark EDM is proportional to 7^ = a.rg(Ad + /i*tan/3) (for the w-quark take tan/5 — > cot/3). 
There are also significant contributions to the neutron EDM from the color dipole moments of the 
quarks, which have the same dependence on the phases at the quark EDMs, and are of the same 
order as the contribution from the induced 9 term described above. The neutron EDM written in 
terms of the quark EDMs and cEDMs (which can be found in lij) is given in pi| . Lastly, there 



is also a small contribution from the three-gluon operator Oq = —^f abc G a GbG c ||22|| . Two-loop 



contributions to the EDMs may be significant at very large values of tan/5 |[23|| . For the values of 
tan j3 we consider, the effect of higher order contributions will be to very slightly shift the allowed 
regions of parameter space. 

For the 199 Hg atom, the electric screening of the electric dipole moments of the atom's con- 
stituents is violated by the finite size of the nucleus and can be conveniently expressed by the 
Schiff moment S, which parameterizes the effective CP-odd interaction between the electron and 
nucleus of spin /. The Schiff moment, in turn, can be induced either through the Schiff moment 
of the valence nucleons or through the breaking of time reversal invariance in the nucleon-nucleon 
interaction, the latter being enhanced by the collective effects in the nucleus, and hence is dom- 
inant. The largest contributions to the Schiff moment, and hence to the EDM of mercury, are 
through the color EDMs of the quarks, and the expression for o?Hg in terms of the quark cEDMs 
can be found in [13]. We emphasize that the mercury, electron and neutron EDMs all depend 
on different combinations of phases, and so simultaneously imposing all three EDM constraints 
excludes a much greater portion of the CP-violating SUSY parameter space than from imposing 
any single constraint. 

Again, the EDMs generated by the SUSY phases are sufficiently small if either 1) the phases 
are very small (< 10~ 2 — 10~ 3 ), or 2) the SUSY masses are very large (O (a few TeV)), or 3) there 
are large cancellations between different contributions to the EDMs, or by a combination of these 
effects. It is condition 3), large cancellations between different contributions to the electric dipole 
moments, which has spurred the greatest interest recently, since it ostensibly allows for both large 
amounts of CP violation and a spectrum which is phenomenologically relevant. Such cancellations 
are not easy to achieve, however. If the phases are 0(1) and the SUSY masses are in the 100 GeV 
range, then the EDM limits will be violated by several orders of magnitude unless very delicate 
cancellations between the various contributions exist. Further, the parameters must be tuned so 
that similarly delicate cancellations occur for electron and neutron and Hg EDMs. 

In some models it is impossible to achieve sufficient cancellations in all three EDMs to permit 
some phases to be 0(1), but even so, the effect of cancellations may still be significant. For 
example, cancellations do occur naturally in mSUGRA models over significant regions of parameter 
space |IJ, 0, [Tj| including in the body of the cosmologically allowed region with my 2 = 
0(100 — 400 GeV). With a .B-type neutralino LSP, large sparticle masses cannot be invoked to 
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permit large phases, due to limits on the LSP relic abundance [|19|, [I], [I8[ , and so the phases, 
in particular 9^, are severely limited by the EDMs. The presence of cancellations relaxes the 
constraints on the phases, but the limit on 9^ remains small, 9^ < 7r/10, unless the sfermions 
(including selectrons) are heavier than O (1 TeV). This upper bound on 9^ is over an order of 
magnitude less restrictive than what would find in the absence of any cancellations, however, and 
thus the effect of cancellations is quite significant. We will discuss the role of cancellations in the 
MSSM in more detail in Section |3.2j . 

Minimal SUGRA is a particularly restrictive model in that there are only two new physical 
phases. In more general models there are more phases, and studies have been made to examine the 
new cancellations which the presence of additional phases allow, particularly in models without 
gaugino mass unification, where the gaugino masses Mj can have independent phases (see e.g. 
H |24|]). The possibility of having both 0(1) phases and reasonable sparticle masses in these models 
has inspired a remarkable number of recent papers exploring the consequences (phenomenological 
and cosmological) of new large sources of CP violation in SUSY models. Most of these analyses 
take the possibility of cancellations as carte blanche to consider all sets of masses and phases, 
without actually imposing the rigorous constraints on the SUSY parameters from the electric 
dipole moments. In the next section we will examine the size of the phases one may reasonably 
expect to satisfy the EDM constraints, and for what sparticle masses, in both mSUGRA and in 
two more general models without gaugino mass unification. In particular, we will study the level 
of tuning required to obtain large phases. 

3 EDM Analysis 

We begin by studying the constraints imposed by the EDMs in three different models: one 
mSUGRA-inspired model with two physical phases, and two models without gaugino mass uni- 
fication, which have seven independent phases each and 15(23) total free parameters, including 
masses. We have done large Monte Carlo studies for each case, evaluating roughly 800,000 pa- 
rameter sets in mSUGRA and 300 and 600 million sets in each of the two more general models 
respectively, and studied the configurations which satisfy the experimental limits on the EDMs. 
Table p] displays the number of parameter sets studied and total number of points satisfying the 
EDM constraints ([1])-®. 

3.1 mSUGRA-inspired Model 

After performing the field redefinitions described in Section |2|, the mSUGRA-like model is specified 
by 6 parameters: three masses (mo,my 2 , and A ), two phases (9^(M X ) and 9a(M x )), and tan/3. 
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Model 


Points run 


Points satisfying EDMs 


fraction satisfying EDMs 


mSUGRA 


8 x 10 5 


10,000 


1.2 x lO" 2 


15 parameter MSSM 


3 x 10 8 


23,650 


7.1 x 10~ 5 


23 parameter MSSM 


6 x 10 8 


13,700 


2.1 x 10" 5 



Table 1: Monte Carlo studies of electric dipole moments in three models. 



Throughout we take m top = 175 GeV. In mSUGRA, once the gaugino, soft scalar masses, A and B- 
terms and phases are given at Mx, they can be evolved using the renormalization group equations 
(RGE) to the electroweak scale. As in common practice, we use the one-loop RGEs for the masses 



and two-loop RGEs for the gauge and Yukawa couplings |p5| . The structure of the equations for 
the gauge couplings, gaugino masses and the diagonal elements of the sfermion masses are such 
that they are entirely real. The evolutions of the however, are more complicated, as the A t 
pick up both real and imaginary contributions. For example, the evolution of A t is given by 

dAt _ 1 ( 16 „2 n/ r ,„2u 13 „2 xjr i 1,2 A ,«j,2 



Sg 2 2 -g( + K + 3hi + 3h 2 t ), (5) 



(Jf s _, v — g* M 3 - 3gi M 2 - - g\ M x + h 2 b A b + 6hfA t j . (4) 

Thus, A t receives real contributions proportional to the gaugino masses Mj and complex contribu- 
tions from the heavy generation Aj, multiplied with the respective Yukawa coupling hi. Since the 
coefficients of the Mi are flavor dependent and the coefficients of the tijAf terms are generation 
dependent, the phases (and magnitudes) of the Ai must therefore be run separately. At one loop, 
the evolution equation for /i is given by 

dfi fi 
~dt ~ 16tt 2 

and the phase of /i does not run. Finally, the B parameter evolves as 

dB 1 
~dl ~ 8^2 

After evolving the parameters to the weak scale, the phase of Higgs superfield H 2 (which gives 
mass to up-type fermions) can be rotated in such a way as to ensure real expectation values for 
the Higgs scalars. The rotation changes the phase of H 2 by an amount —{0^ + 8b)- Not only is 
the phase of fi now fixed at 9^ = —6b, but also the initial phase of fi is physically irrelevant as 
it is canceled by the rotation. As emphasized in [26, 27], a large phase in A will induce a phase 
in B (Eq. (^)), and hence in fx , after the vevs are made real. Therefore even if fi and B are 
both real at Mx, if Qa is large, the value of Q^iMz) contributing to the EDMs may be large (of 
course this is completely equivalent to keeping the Higgs vevs complex and fi real). Since B is a 
free parameter, which is typically determined by the conditions of correct electroweak symmetry 



-3gl M 2 - gf M l + h T A T + 3h\A h + 3h 2 t A t . (6) 
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breaking, B(Mx) can be chosen so that Bfi(Mz) is nearly real, yielding a small 9 VL (Mz) after the 
Higgs rotation; however, this can involve a significant fine-tuning if 9 a is large f2"7}| . In practice, 
this tuning is typically not worse than at the level of 10% |Tj|, but it must be emphasized that 
this is a tuning over- and- above the tuning discussed below. Alternatively, if 9b{Mx) is taken real 
2T)[ and 9a is large, or if 9b{Mx) is arbitrary, then the additional fine tuning lies in adjusting the 



original 9® close to —9b(Mz)- 

To study the phases permitted by the current experimental limits on the neutron, 

electron and 199 Hg EDMs, we perform a Monte Carlo studies, sampling the 6 model parameters 
in the following ranges: 

< 9^9 A < 2tt 
50 GeV < m V2 < 1 TeV 
GeV < A < 1 TeV 
GeV < m < 1 TeV 

2 < tan/3 < 10 

and keeping parameter sets satisfying (0)-©. Limits from particle searches have not been imposed. 
However, the lower limit on my 2 was chosen to remove most of the area excluded by the current 
LEP chargino mass limit of 103 GeV ||28|| . Most of the solutions with tan/3 < 3 are excluded by 



the Higgs searches at LEP||28||. Cosmological limits on the relic abundance of the LSP neutralino, 
which tend to exclude regions with mo ^ 200 GeV, have also not been included, we will comment 
on these regions later. 

The results for our Monte Carlo scan of 800,000 mSUGRA parameter sets are summarized in 

Fig. [| where we display the 10,000 setsQ satisfying the three EDM constraints (|1])-(|3|). Fig. |3]a 

shows the allowed configurations of {9^, 9a}- The first obvious point to note is that the range of 

9/j, is severely limited, with no events in this sample having > 0.17T (mod tt), while the phase 

9a can take any value. This is due to the fact that the chargino contribution to the electron EDM 

is typically dominant and depends only on 9^; therefore for large the neutralino exchange 

piece cannot provide sufficient cancellation, regardless of the values of 9a- In Figs. |||b-|3|d, we 

display the value of 9^ versus mi/ 2 , m and tan/3, respectively. We see that there are many more 

solutions for large > 500 GeV, where the heavy spectrum reduces the individual contributions 

to the EDMs. We note, however, that almost all such parameter sets lead to a neutralino relic 

abundance &zh 2 > 0.3 [ 15fl , implying a universe younger than 12 billion years, in contradiction 

to observational evidence . In Fig. ^d, we see that most of the large 9^ solutions occur also for 

tan/3 < 3, which in mSUGRA yield a Higgs scalar too light to be compatible with the negative 

results from Higgs searches at LEP [p8 |. Indeed, if we consider only those parameter sets with 

mi/2 < 500 GeV, mo < 200 GeV and tan (3 > 3, we find no solutions with 9^ > 7r/20. 

2 To reduce the file size in this and following plots, we do not display points lying underneath covered regions, 
particularly at low 8^. This produces no visible effect on the plots. 
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3.2 15-parameter MSSM 

We next consider a more general model which has additional independent phases, and which 
therefore has greater opportunity for cancellations. We no longer require gaugino mass unification. 
Therefore, one of the three gaugino masses may still be made real by an R-rotation, but the other 
two, which we take to be Mi and M 3 , may be complex, with phases 0i and 3 respectively. 
We additionally allow independent phases in Ad, A u , A e and A t , and along with 6^, this gives 7 
independent phases. The other 8 parameters are tan/3, plus the masses: |Mj|, % — 1 ... 3, a common 
\Aj\ = \A\,i = e,d,u,t, and the sfermion mass parameters mf ,m| . We fix the remaining 
masses by the approximate relations m~ = m\ , m| = m| + 0.6M|, and = m? + 0.5Mf. 
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Figure 3: mSUGRA solutions satisfying the experimental limits on the electron, neutron and 199 H^ 
EDMs. 
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We then sample the regions 
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and look for parameter sets satisfying (|T])-(^|). We cut out the (few) remaining sets with m x ± < 
103 GeV, and we flag the sets with m h < 113 GeV, as discussed below. 

Again, because of the additional possibilities for cancellations in this parameterization, we 
expect solutions with larger values of the phases which still satisfy the EDM constraints (|l|)-(|3|). 
This is borne out in Figs. |J-|6|, where we display the results for our Monte Carlo scan of roughly 
300 million points in the 15-parameter MSSM. In the top two panels of Fig. |] we display the 
23,500 sets allowed by the EDMs, in the {9^,(pi}, and {9^,(fi 3 } planes, respectively. Although by 
far the greatest events of solutions have 9^ < vr/10, there is now a visible swath that extends 
to large values of 9^. The range of </>i and </>3 are unconstrained, as are the ranges of the 9^ 
(not displayed). The presence of the large 9^ solutions relies on having large sparticle masses, 
significant cancellations between contributions to the EDMs, or small tan j3, and very typically a 
combination of the above. 

In Fig. |] we display 9^ versus tan (3 for the parameter sets, and we see clearly the weakened 
constraint on 9^ at low tan f3. However, most of the low tan (3 solutions are actually experimentally 
excluded by the unsuccessful Higgs searches at LEP2. To be conservative, we take the stop soft 
masses to be independent parameters and set them to 1 TeV, take = 1 TeV, and we compute 
the light Higgs mass as a function of tan/? in the maximal mixing scenario using Ref. |[29|| . The 
light (green) points in all four panels of Fig. [| yield < 113 GeV (corresponding to tan/? < 2.4). 
Hence, they would require large phases and further tuning in the Higgs sector not to be excluded 
by LEP f28fl . In Fig. [5], we display 9^ as a function of the light {xf) and heavy (xt) chargino 
masses. The chargino exchange contribution to the fermion EDMs and cEDMs fl"4] , ||| is suppressed 
both for highly split charginos and for a heavy chargino spectrum. We see that if neither chargino 
is heavy {m x ± < 500 GeV), 9^ is strongly constrained, whereas if at least one of the charginos is 
heavy, 9^ can potentially be large. In particular, having a spectrum with only one light chargino 
does not by itself forbid 9^ > n/10. 

The significance of cancellations on the allowed range of 9 ^ is seen in Fig. |[ In the first panel 
we display the value of 9^ versus the ratio the the chargino (d^ ) and neutralino {d^) exchange 
contributions to the electron EDM. We see the clustering of the large 9^ solutions near to the 
cancellation point d^/d^ = —1. The vertical lines correspond to an electron EDM a factor of 5 
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Figure 4: Parameter sets in the 15-parameter MSSM satisfying the experimental limits on the 
electron, neutron and 199 Hg EDMs. Open circles suffer from parameter tuning AX/X worse 
than 1% (see the text). Light (green) dots correspond to configurations with a light Higgs m h < 
113 GeV. The tuning parameter AX/X is defined in the text, and corresponds to the maximum 
variation that the point survives. 
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smaller than the larger of d^ and , due to cancellations. This demonstrates that the sparticles 
must be much heavier than those allowed by our range of parameters (> 1 TeV) in order for the 
EDM constraints to be satisfied for large 9^ in the absence of any cancellations. In practice, the 
large 9^ solutions have both some sparticle masses near the upper end of their ranges and a large 
degree of cancellations. In the second panel of Fig. || we display 9^ as a function of the sfermion 
mass m e - R . Perhaps counter- intuitively, the large 9 ^ solutions tend to lie at lower rrie R , rather than 
near the upper end of their range. However, this is again due to the necessity for cancellations: 
mi ~ 1 TeV is not large enough to sufficiently turn off the SUSY contributions to the electron 
EDM for large 9^, and the necessary cancellations only occur for smaller slepton masses. 

One goal of this paper is to study the extent to which the large phase solutions require a 
significant tuning of the model parameters in order to fall into the regions of EDM cancellations. 
The tuning measure we employ is simple but intuitive. For every parameter set we find which 
satisfies the EDM constraints, we perform the following analysis. We begin by varying all the 
input parameters (one at a time) by AX/X = ±0.5% and see if the EDM limits are still satisfied 
for all the test parameter sets. If they are, we then try varying all the input parameters by ±1%, 
and so on, until we find the smallest percentage change for which the configuration violates one 
of the EDM bounds. This gives a sense of the local "size" of the allowed parameter region. The 
results for the 15-parameter MSSM scan are displayed in the final panel of Fig. [|, where we have 
stepped through 0.5%, 1%, 2%, 5%, 10%, 20% and 50% changes, modulo tc in the case of phases, 
and the last percentage change to the underlying parameters that successfully satisfies the EDM 




Figure 5: For the 15-parameter MSSM scan, the dependence of the solutions on the light and 
heavy chargino masses. Open circles suffer from parameter tuning AX/X worse than 1% (see the 
text). Light (green) dots correspond to configurations with a light Higgs < 113GeV. 



11 



bounds is plotted versus 9^. Points which do not survive a 1% change are not plotted in the last 
panel of Fig. |], although they are plotted in all the other figures. It is clear that the large # M 
solutions require more tuning that the low 9^ solutions. Overall, we find that {27%, 46%, 74%} of 
the (large plus small 9^) solutions require worse than {1%, 2%, 5%} tuning (corresponding to {the 
unplotted points, points at 0.01, points at 0.02} in the last panel of Fig. |]). The tuning level of 
the solutions can also be seen in the other panels of Figs. [|, where the open circles have tuning 
worse than at the level of 1%. 

We emphasize the importance of having a third independent constraint on the phases, from the 
199 Hg EDM. Fig. |7] shows the result of a smaller scan of the 15-parameter space in which we only 
impose the limits from the electron and neutron EDMs and plot 9^ vs. tan/3 for the resulting set. 
The light (green) points are those which would be forbidden by the further imposition of the 199 Hg 
EDM constraints, while the dark (black) points are those which satisfy the 199 Hg EDM constraints. 
The open circles yield rrih < 113 GeV. We find that of the 174,000 points satisfying the neutron 
and electron EDM constraints (of which only 50,000 are plotted), only 4700, or 2.7%, additionally 
satisfy the 199 Hg EDM constraint. The remaining points at large 9 ^ typically have a very heavy 
squark (> 800 GeV) or two heavy charginos (> 500 GeV). Clearly, ignoring |3(| the 199 Hg EDM 
constraints allows many configurations that are experimentally forbidden, particularly for low to 
moderate values of the masses, and for larger values of tan (3 combined with large phases, where 
the regions satisfying the electron and neutron EDM constraints are small. Given the significant 
effect of the mercury constraint, an improved calculation of the strength of the T-odd nuclear 




Figure 6: For the 15-parameter MSSM scan, the effect of cancellations on the allowed range of 
9^. Open circles suffer from parameter tuning AX/X worse than 1% (see the text). Light (green) 
dots correspond to configurations with a light Higgs < 113 GeV. 
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forces |5T| and Schiff moment of the mercury nucleus [E2J will likewise be important. 

We also emphasize the significance of the recent improved calculation of the neutron EDM 



21]. In Fig. |8| we display the values of the neutron EDM as computed using QCD sum rules in 
|2"T] against those estimated using naive dimensional analysis (NDA), as in Here we have 
plotted just those sets for which the neutron EDM constraint is satisfied according to either one 
or both of the calculational methods. We see that NDA typically overestimates the neutron EDM 
by roughly a factor of two, although for solutions near regions of cancellations, the discrepancy 
between the two can be much greater. 



3.3 23-parameter MSSM 

We lastly consider a very general model with 23 free parameters. We now allow the left and right 
sfermion masses to vary independently, and we take independent stop masses in computing the 
neutron EDM. We also take the phase and magnitude of all the A^i = e,d,u,t as independent, 
giving a total of 7 phases, 15 mass parameters and tan (3 as the free parameters of the model. We 
perform a Monte Carlo scan of roughly 600 million parameter sets, over the following ranges: 




0.5 1 1.5 2 2 4 6 8 10 

fy^n tan (3 



Figure 7: A scan of the 15-parameter MSSM imposing only the electron and neutron EDM con- 
straints. Light (green) points are forbidden by the 199 Hg EDM constraints, while the dark (black) 
points satisfy the 199 Hg EDM constraints. The open circles have < 113 GeV. 
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d N1 in units of 10" e cm 




Figure 8: Computed values for the neutron EDM, using QCD sum rules (y-axis) and naive dimen- 
sional analysis (x-axis), in units of 10~ 25 e • cm. 
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In Fig. | we display the 10,000 parameter sets satisfying the EDM constraints (TJ)-(0)- The results 
are very similar to those in the case of the 15-parameter MSSM. 



4 Mass and Cross Section Measurements 



We now turn to the actual determination of the phase parameters at a future Linear Collider [[33 
In an earlier work [|[ it was shown that large phase values, as well as the real parameters in the 
neutralino/chargino sector, can be extracted to high accuracy from measured masses and cross 
sections, using a global fit. Their extraction from a much reduced set of observables, such as the 
neutralino and chargino masses alone, appears to be impossible, due to experimental uncertainties 
propagating from the measurement of masses and cross sections into the fitted phase parameters. 
Only if the uncertainty is considerably smaller than the actual phase value, do we regard the phase 
as being observable. 

4.1 Complete Set of Observables 

The central scenario in our analysis is a generic set of MSSM parameters derived from the unified 
model parameters m%/2 = 200 GeV, m = 100 GeV, A Q = 0, tan/3 = 4, and ji > 0. The 
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Figure 9: Parameter sets in the 23-parameter MSSM satisfying the experimental limits on the 
electron, neutron and 199 Hg EDMs. Open circles suffer from parameter tuning AX/X worse 
than 1% (see the text). Light (green) dots correspond to configurations with a light Higgs m h < 
113 GeV. The tuning parameter AX/X is defined in the text, and corresponds to the maximum 
variation that the point survives.. 
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neutralino/chargino mass parameters are \M\\ = 83 GeV, M2 = 165 GeV and /x = 310 GeV, and 
the corresponding slepton masses are rrie L = 180 GeV, m e - fl = 132 GeV and mo = 166 GeV. The 
masses of the neutralinos and charginos are given in Table |2|. For the phase values <pi = 9^ = n/10, 
this scenario was investigated in our previous analysis: the error on the fitted mass parameters 
|Mi|, M 2 , \fi\ is smaller than 1 GeV, and the error on the extraction of tan j3 is less than 10%. The 
RMS of the phase determination depends on the energy and the luminosity of the collider. For 
a (500 GeV, 500 fbT 1 ) machine we obtain (fii/n = 0.1 ± 0.03 and = 0.1 ± 0.05, whereas for 
(ITeV, lOOOfb" 1 ), the errors are fa/ir = 0.1 ±0.05 and 0Jir = 0.1 ±0.060. The t-channel slepton 
masses are assumed to be measured in threshold scans, and the propagation of the uncertainty 
of their mass measurements into the phases can be neglected. We find that the statistical errors 
on the phase determination A0i ~ 0.037T and A6* M ~ 0.057T are essentially independent of the 
central values of the phases. This implies that even with a maximal set of observables, phase 
values smaller than 7r/10 are hidden by experimental errors and are therefore unobservable: CP 
phases of that size could, from a collider phenomenology point of view, as well be zero. 

A naive view of the inversion problem is given by the graphs in Fig. [H| We fix the parameters in 
the neutralino and chargino mass matrices and plot the production cross section a(e + e~ — > X1X2) 
and the branching fraction BR (%!] — > ^e + e~) as a function of the phase <fr\. The dependence on 
<pi is strong in both observables and an analysis could be straightforward. However, the plots also 
show that the physical masses in the process are not constant with varying phases, either. The 
main problem of the analysis therefore becomes the separation of direct effects of the phase and 
indirect effects, where the phase changes the physical neutralino masses and these affect the cross 
sections and branching fractions. 

3 For both colliders the masses are determined in a threshold scan. For higher energy and fixed masses the cross 
sections decrease. Increasing the luminosity does not compensate for the increase in energy in this case. 
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Hence, a key ingredient of the full analysis 0] as well as of the one presented in this paper is the 



measurement of masses through threshold scans |33|. We use the error estimates for an mSUGRA 



scenario with m = 100 GeV, m 1/2 = 200 GeV and tan/3 = 4. They are 0.05/0.07/0.3/0.6 GeV for 
the neutralino masses and 0.035/0.25 GeV for the chargino masses ]3"4| . We furthermore assume 
the uncertainty of cross section measurements to be purely statistical. This limits the pull of small 
cross sections. The invisible cross section for the production of two lightest neutralinos (LSPs) is 
not part of our sample. For a high energy collider, the full analysis includes twelve cross section 
and six mass observables, determining six model parameters: \Mi\, M 2 , \fi\, tan/5, <pi, 

4.2 Inclusion of the EDM limits 

To cover a large range of slepton masses, we consider two modifications of the central scenario: 
(1) the first generation sleptons are still light but just escape detection at a low energy collider. 
The so increased slepton masses correspond to having their mo = 200 GeV, m 1( / 2 = 200 GeV, i.e. 
nie L = 250 GeV, = 218 GeV and ma = 240 GeV. (2) All sleptons decouple from the theory 
entirely, which corresponds to masses above 0(1 TeV). The cross sections accessible at a 400 GeV 
collider are given in Table ^| for all three sets of lepton masses. This collider energy enables 
us to observe one higgsino type neutralino directly, whereas the other two higgsino states are 
kinematically inaccessible. The cross sections are dominated by the t-channel process due to the 
gaugino nature of the light neutralino. In contrast, the production of two light charginos involves 
both diagrams: for light sleptons, the t-channel graph is large whereas in the decoupling limit the 
s-channel contribution dominates; for intermediate slepton masses destructive interference reduces 
the cross section by almost a factor of two. Still, the chargino cross section is by far the largest. 
If the lightest chargino mass is known through a threshold scan, this signature should serve to 





[GeV] 




a tot [fb], i light 


Ctot [ft"]' 1 


? heavy 


ctot [ft^ 1 


I decoupled 


X? 

xi 

xl 
xl 
xt 
xt 


77.3 
142.6 
315.5 
343.0 
140.7 
431.6 


(xSx?) 

xtixl 

xlxl 

xixi 

xtxi 


320.0 
87.8 
6.4 
107.8 
806.9 


186.6 
55.3 
3.8 
64.5 

570.8 


0.09 
0.28 
1.7 
0.16 
854.6 



Table 2: Neutralino/chargino masses and cross sections at a 400 TeV Linear Collider for the 
three scenarios under consideration. The production cross section of two LSPs is not part of the 
set of observables. 
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determine the mass of the sneutrino with high precision |35| . 

Using the results of the EDM analysis allows us to modify the full analysis described in the 
previous section: we know that the phase 9^ has to be smaller than vr/10, after taking into account 
three conditions: (1) the experimental limits on the electron, neutron and mercury EDM have to 
be respected simultaneously; (2) the degree of fine tuning is limited by requiring that the solutions 
be stable with respect to changes of AX/X = 1% in all model parameters; (3) a minimal set of 
final states is produced at a Linear Collider with a fixed design energy, e.g. 400 GeV. We will 
show that this minimal set of observables has to include a handle on one higgsino component, to 
allow for a determination of fi. Under these assumptions, we try to determine CP-violating phases 
at a comparably low energy collider, which would not cover all neutralino/chargino and slepton 
thresholds. The reduced set of observables is supplemented by EDM constraints: tan (3 has to be 
small (< 10), and may be known from the Higgs sector |3(| if the Higgs bosons are sufficiently light. 
Since 9^ turns out to be constrained to be smaller than its minimum visible value 7r/10 0, we fix 
it to zero. The phase <pi can be treated as independent of all other parameters in the neutralino 
and chargino mass matrices, since parameters correlated to 0i by the EDM constraints, like 03 or 
4>Ai do not appear in this sector. 

To determine the phases we modify our full analysis: the set of observables is first reduced 
to the three lightest neutralinos and the lightest chargino. Producing them in pairs yields four 
cross sections, i.e. eight independent observables. In the fit, we determine the real parameters 
| Mi |, M2, /i, 4>\ as well as two t-channel selectron masses. For a low energy collider, the selectron 
masses cannot be expected to be measured in threshold scans, and this leaves us with altogether 
six unknown model parameters. Using SU(2) symmetry the sneutrino mass can be related to the 
left handed selectron mass: 

m~ = m\ L + cos 2 9 W cos 2(3 . (7) 

In contrast to the complete parameter fit, the extraction is now limited by the number of fitted 
parameters, and the \ 2 distribution of the best fit might be flat in certain parameters. This 
makes it technically difficult to add tan (3 to the set of fitted parameters, and we have to rely on a 
measurement in the Higgs sector. However, in principle one can extend the number of observables 
by non-trivial distributions, asymmetries or additional cross sections, and one can choose a fitting 
algorithm better suited for the problem to include tan f3 in the fit || . 

4.3 Statistical Uncertainties 

The inclusion of experimental errors follows the same path as in the previous analysis |4|]: we 
assume Gaussian errors for the measured masses and cross sections and define smeared pseudo- 
measurements. First we define a set of 'true' model parameters \Mi\, M 2 ... They predict a set 
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of 'true' observables (masses and cross sections), all of which are assumed to have a Gaussian 
error distribution with a known width |[34| . Using these distributions, we randomly vary 10000 
pseudo-measurements of the set of observables. These sets become slightly inconsistent, but using 
a global fit we can extract the central value of every model parameter and obtain a distribution 
for each of themQ. In the fit we minimize x 2 — J2i (^reconstr,i ~ ^meas.i) 2 / 'ef , where x reconstr are the 
reconstructed observables and x meas are the smeared pseudo-measurements. The error on the cross 
sections is given as a function of the luminosity by e; = cr;/(e£), where the efficiency is assumed 
to be e = 10%. If the central value of the distribution for a given parameter agrees with its 
'true' value, the statistical treatment is justified, and the width of the distribution describes the 
migration of observational errors into the model parameters. The final curve does not necessarily 
have to be Gaussian, since correlations, together with the range of the fit, can alter the shape of 
the curve. Hence, we quote the RMS value instead of the standard deviation of a fitted Gaussian 
distribution. We note that for some scenarios there exist several x 2 minima in different parameter 
regions. This is a technical complication, and the minima should be distinguishable through the 
actual values of x 2 ■ 111 these cases, we limit our range of fitted parameters, making sure that the 
range is much bigger than the distribution we finally obtain. 



The distribution of best fits to a set of 10000 pseudo-measurements is given in Fig. [11]. In 
contrast to the earlier analysis we assume 9^ = and rely on a known value of tan/3. The 
sleptons are kinematically inaccessible, and their masses have to be indirectly determined from 
the cross section measurement. Comparing the errors on the three mass parameters \Mi\,M2,fi 
we notice the striking accuracy of the measurement of M 2 . This reflects the small number of 
unknown parameters in the chargino mass matrix and in the xfxl production cross section: 
the phase 9^ and tan/3 are both fixed, and only M 2 and \i and the sneutrino mass have to be 
determined from the chargino mass and cross section. However, if one is able to determine the 
cross sections for polarized electrons and positrons and separate the s and t-channel contributions, 
this accuracy will improve even further. On the other hand, if one is forced to determine tan/? 
from the neutralino/chargino sector, the error on M 2 will be 0(1 GeV) again. In contrast to 
the full analysis, the fitted values of \M\\ and are now correlated. Correlations like that are 
a general feature of a smaller set of observables and a smaller dimension of the fit. The error 
on the determination of the slepton masses is not symmetric: the left selectron mass is tied to 
the sneutrino mass and therefore extracted from the well determined chargino sector. The right 
selectron mass, in contrast, has to be determined from the neutral gaugino cross sections, which 
are smaller by almost one order of magnitude. 

The distribution of best fits for <pi exhibits a similar shape to the result of the full analysis. For 



4 Extracting these values from the low energy parameters without the smearing would merely serve as a check 
of the fitting program. 
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large luminosity of 500 fbT 1 it approaches a Gaussian with a RMS of 0.032 and a fitted standard 
deviation of 0.028. At smaller luminosity, this distribution broadens, e.g. to a RMS value of 0.037 
for 300fb _1 , as given in Table |3|. From the cross sections in Table |2] one can see that there is 
a minimal luminosity at which the number of events in (X1X3) becomes marginal. In this case 
there will be hardly any information on \i left in the sample. The same problem arises once the 
energy is too small to produce a state with significant higgsino content. If the higgsinos are almost 
entirely decoupled, the uncertainty on the real mass parameters becomes too large to extract the 
phases. Technically, the correct minimum in y 2 will not be found then by the fitting procedure. In 



Figure [T0| we see that for many observables the point 0i = is extremal, i. e. a considerable number 



of best fits will give this result. This occurs for example in the case of low luminosity, Fig. [IT]. On 
the other hand we notice that in the set of observables under consideration, the gaugino sector 
is still over-determined: except for the less important slepton masses the X2X2 production cross 
section carries similar information as the much larger XiXi cross section. Removing the X2X2 
cross section from the sample increases the RMS value given in Table |] to 0.036 ; however, in 
most models X2X2 wu ^ De visible at any collider that produces pairs of light charginos. 



83.5 



(0.4TeV, 500 fb" 1 ) (0.4TeV, 300 fb" 1 ) 
M,[GeV] 



83 



82.5 




o.i 



0.2 



311 



310 



309 



. ^[GeV] 


















' M 2 [GeV] 

If , 



164 



165 



166 



260 - 



240 



M(e L )[GeV] 






BBSS©."' . 




• .v- r * •*• . 


, , , ft, 


• M(e R )[GeV] 


200 


250 



10000 best fits 



2000 




(|) 1 /7t=0.098±0.032 

(11/71=0.098+0.037 



0.2 



Figure 11: Parameters obtained from best fits to 10000 pseudo- measurements. The central scenario 
is indicated by arrows on the axes. Only one Higgsino X3 is assumed to be observed. The black 
and grey (green) points correspond to two different collider designs. 
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4.4 Systematic Uncertainties 

As mentioned above the errors derived above do not take into account any systematical errors 
from using a wrong value e.g. of tan/3. Neither do they include an estimate what would happen 
if # M were actually small but non-zero. Both effects would lead to systematic errors and thereby 
to wrong central values of the 10000 fits. In Table § we show that fitting a set of parameters 
with a 'true' value # M = vr/20, but assuming 6^ = hardly affects the fit. This observation serves 
as a consistency check of the ansatz, namely that a small 9^ has no effect on the extraction of 
<pi and certainly not on the dominating real parameters in the mass matrices. Analyzing the set 
of observables with a wrongly assumed values of tan/3 = 3.8 or tan/3 = 4.2 instead of a 'true' 
value of 4.0 leads to the systematic error shown in Fig [12|: the spread of the best fits in the mass 
parameters as well as in the phase is very similar to the case where both parameters, # M and tan f3 
are input correctly, but the central values are shifted. The central values for wrongly assumed 
values tan/3 = 3.8 and 4.2 become fa/ir = 0.125 ± 0.033 and fa/ it = 0.082 ± 0.029, i.e. they are 
shifted by half the RMS value. 

To estimate the dependence of the analysis on the details of the underlying model, we decouple 
the sleptons from all cross sections. From the cross sections in Table 0, we see that the real mass 




Figure 12: The effect of wrongly assumed values for tan/5 on the extraction of the model param- 
eters. The 'true' value for tan (3 is 4. 
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parameters dominated by the physical chargino/neutralino masses and the chargino cross section 
should be extracted as precisely as in the case of light sleptons. However, the determination of 
the phase </>i relies on the measurement of the X1X2 an d the X1X3 cross section, where the Z boson 
in the s-channel production process for neutralinos only couples to the higgsino fraction of the 
physical states. Since in our central scenario the light neutralinos are mainly gauginos, both of 
these cross sections decrease with increasing slepton masses. This leads to a RMS value of 0.044 
of the 0i distribution, as seen in Table |[ 

An important feature of the complete analysis j| is that the error on the <pi measurement is 
independent of the actual value of <\>\. This allows us to derive lower limits on the size of observable 
phases, which are 0i > 7r/10 and # M > 7r/10 for the complete analysis. As presented in Table H, 
we obtain RMS=0.030 for a central value of <p\ — tt/2. In Figure [Hj the cross section, as well as 
the branching fraction, tends to vary strongly around phase values of tt/2, while they become flat 
and extremal for positive or negative real values of Mi. This is reflected in the slightly smaller 
RMS value for </>i = tt/2 as compared to 0i = 7r/10. 

From the cross sections in Table |2] it is obvious that the determination of </>i from the given set 
of cross sections improves, once fi takes a value closer to the gaugino mass parameters. Reducing 
it by a factor of two to fi = 150 GeV leads to a much lighter higgsino with m^o = 158 GeV and an 
increased cross section for X1X3 production. The changes in the masses of the gaugino dominated 
light neutralinos are small, but the error on the </>i measurement improves to 0.014. 

Although it is hard to accommodate with the EDM constraints, we investigate a scenario with 
large tan f3 = 30: from the full analysis we expect the determination of tan f3 to be less precise than 
for smaller values. Since all production cross sections are smaller by up to 50%, the extraction 
of the slepton masses becomes increasingly difficult. From the analysis for small values of tan/3, 
however, we know that the error on the determination if the right selectron mass is larger than 
the mass difference between the two selectrons. One way of improving the fit for large tan/3 is 
therefore to assume that the selectrons are mass degenerate. With this assumption we are able to 
determine the phase </>i = 0.11 ± 0.03. 

4.5 Minimal Set of Observables 

Since the EDM constraints on ^ < 7r/10 require the chargino sector to be essentially CP- 
conserving, we can define a minimal set of observables sufficient to extract the phases: in the 
CP-conserving case one could be able to extract sufficiently precisely from a set of the mass, 
cross section and asymmetries in XiXi production analytically || [37|. We use this result to 
investigate the phases in the neutral gaugino sector, including the masses and production cross 
sections of Xi an d X®- m our central scenario (Table ^) this requires a minimum collider energy 
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of 320 GeV. To the set of three cross sections and three masses we fit the neutralino parameters 
|Mi|, M2, 0i, for given \fi\, and obtain an error on the determination of <fi\ of 0.027, i.e. similar 
to case including one higgsino. Adding tan (3 to the fitted parameters increases the error on <pi to 
RMS= 0.032, as can be seen in Fig. [13]. As shown before, the inclusion of slepton masses would not 
change this result significantly. It only requires a larger set of observables. A wrong measurement 
of |/i|, however, would lead to a systematic error of the phase determination. We estimate this by 
assuming a wrong extraction of = 305 GeV and 315 GeV for a 'true' value \fi\ = 310 GeV. In 
Fig. 13 the extracted values of |Mi| and M 2 are shifted systematically, as is the central value of 
the 10000 pseudo-measurements of <\>\. But the shift for a 5 GeV mismeasurement of |/i| is only 
by half the statistical error of the phase measurement. We therefore conclude that it is possible 
to extract the phase (f>i from a minimal set of parameters, if the error on the determination of \n\ 



does not exceed a few percent [37 
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Figure 13: Top: Fit results for a known value of fi = 310 GeV. The grey (green) points correspond 
to tan j3 being added to the fitted parameters. Bottom: Systematic effect from a wrong input of 
\n\ for known tan/?. 
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j 1 / / 

central (pi /it 


RMS 


(1) 


complete analysis (500 GeV, 500 fb" 1 ) [| 


0.097 


0.030 


(2) 


one phase 0i, free sleptons, fixed tan (3 (500 GeV, 500 fb" 1 ) 


0.099 


0.034 


(3) 


as (2) 


, but only one higgsino xl (400 GeV, 500 fb 1 ) 


0.098 


0.032 




as (3) 


, but (300 fb" 1 ) 


0.098 


0.037 
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, but without X2X2 production 


0.098 


0.036 




as (3) 


, but 'true' M = tt/20 


0.100 


0.032 




as (3) 


, but 'true' 9^ = n/20 and wrongly assumed tan/3 = 3.8 


0.125 


0.033 




as (3) 


, but 'true' 9^ = n/20 and wrongly assumed tan/5 = 4.2 


0.082 
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, but sleptons decoupled 
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0.044 
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, but large phase value </>i = n/2 
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as (3) 


, but light higgsinos /1 = 150 GeV (300 GeV, 500 fb" 1 ) 


0.100 


0.014 
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only £ 


^auginos (320 GeV, 300 fb" 1 ) 


0.098 
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, but 500 fb" 1 
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0.021 




as (4) 


, but tan/? fitted 


0.098 


0.032 




as (4) 


, but wrongly assumed \i = 305 GeV 


0.111 


0.024 




as (4) 


, but wrongly assumed \i = 315 GeV 


0.086 
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Table 3: Central and RMS values for the distribution of </>i values, fitted from the 10000 pseudo- 
experiments, with a "true" value of 0i = 0.17T (except where indicated). 



5 CP-Odd Variables and Asymmetries 

To unambiguously identify an effect of CP violation, one needs to construct a "CP-odd variable" , 
whose expectation value vanishes if CP is conserved. At an e + e~ linear collider, the initial state 
can be made a CP eigenstate, given the CP transformation relation 

e~(ai,p) e + (a 2 , -p) e~(-a 2 ,p) e + (-cri, -p), (8) 

where <j{ is the fermion helicity. We consider a specific process 

e-(pi)e + ( P2 ) -> xlxl - ^(qi)4(Q2)x° 1 X° 1 . (9) 

Denote a helicity matrix element by A / i (Tl o- 2 (9i, 9*2) where <j\ (02) is the helicity of the initial 
state electron (positron), which coincides with the longitudinal beam polarization; qj (q 2 ) denotes 
the momentum of the final state fermion (anti- fermion). For the process of Eq. (|J), only two 
combinations of the helicity amplitude contribute. It is easy to show that under CP 

transformation, 

M-+(q\, q 2 ) => M-+(-q 2 , -ft), (10) 
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and .M+_ transforms similarly. If CP is conserved in the reaction, relation ( |I~0| ) becomes an 
equality. This argument is applicable for unpolarized or transversely polarized beams as well. 

One can construct CP-odd kinematical variables to test the CP property of the reaction. We 
consider the following three angles defined as 

gt ■ q+ (pi x g_) • (fi x q 2 ) gt ■ (gi x q 2 ) 
cos 9 + = , COSfeL = -— ^ — -, cos 9 e = . - . . _ ^-r, (11) 

IpiI Ipi x q~\ m x 92 1 ' M 1 91 x 92I 

where q + — q\ + q 2 and <f_ = q\ — g*2- It is easy to verify that all the three variables are CP-odd 
under final state CP transformations. We can then construct "forward-backward" asymmetries 

A FB = a F -a B = f 1 -^-dcosO - f° -^-dcosO, (12) 



o dcos9 J-i dcos9 

with respect to a CP-odd angle 9. 

Of the four Feynman diagrams contributing to process (|9|), there are two diagrams that contain 
explicit CP-violating phases, one from s-channel Z exchange and one from a t-channel selectron 
exchange. If \x > \Mi\,M 2 {fi < \M\\,M 2 ), then in most part of the parameter space the con- 
tribution from the selectron (Z) exchange diagram is dominant, thus makes the CP asymmetry 
from the interference small. We scanned the parameter space in (0 1; 6^), the asymmetry obtained 
is typically 0.1 — 1.5%. The beam polarizations do not improve the situation significantly. As 
an example, for \M X \ = 80 GeV, M 2 = 200 GeV, n = 275 GeV, tan/3 = 4.0 and m gR = 165 GeV, 
while 0i = 0.907T and 9^ = 0.257T the asymmetry from cos 6^ appears to be about 1.0%, with 
an asymmetry rate of 2.7 fb. The asymmetries for other variables in Eq. fllT|) are comparable. 
Naively, such an asymmetry at a high luminosity linear collider of 200 fb -1 would lead to a 2a 
statistical effect, as estimated in |30[. However, one would have to keep the systematics of the 



asymmetry measurements well below a percent level, in order to establish a positive observation. 
We thus consider it very challenging to experimentally observe this rather small asymmetry. 



6 Conclusions/Outlook 

We have shown that the current experimental limits on the neutron, electron and 199 Hg electric 
dipole moments strongly constrain general SUSY models with several CP-violating phases, even 
in the presence of strong cancellations between the various SUSY contributions to the EDMs. 
Although it is only 9 ^ which is typically constrained to be small, in models in which the phase of 
B is correlated to the other (large) phases in the model, this translates into a tight restriction on 
some combination of (large) phases, which may not be natural in any given model. The next year 
will see significant improvements in the experimental limit on the electron EDM |38| , and work on 
improving the mercury EDM limits continues as well. Further down the road, new experimental 
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techniques, such as using diatomic molecules to measure the electron EDM []55f| , or studies of the 
effects of CP violation on other systems, such as (bb) pro duct ion j40], will probe additional CP 
violation in the MSSM to still more sensitive levels. 

These constraints on CP- violating phases are important for studies at a future Linear Collider. 
It was shown that one can always extract the phases 0i and 9^ from a complete set of masses 
and cross sections for neutralino and chargino pair production [||]. The new EDM constraints 
essentially require 9^ to be too small to be measured at a Linear Collider. In the region where 
charginos and neutralinos are visible at a Linear Collider tan/3 is preferably small, i.e. rendering its 
measurement accessible through Higgs production p6 |. Using the above constraints we extracted 
the finite phase (f>i from a reduced set of masses and cross sections and find that a minimal set 
of observables is limited by the presence of a non- vanishing Higgsino component in the final state 
particles. Since a lower energy (400 GeV) Linear Collider might not be able to produce slepton 
pairs we show that the t-channel masses in chargino and neutralino production can be fitted, 
their uncertainty does not interfere with the <f>\ phase measurement. Pair production of light 
charginos |I] might provide us with precise indirect information on the size of |/i| ||37|| . We showed 
that in this case that we do not need to produce any chargino directly to determine the phase 4>i 
from a truly minimal set of (gaugino) observables. For this set of observables we can even measure 
tan /3 and easily cross check a tan (3 measurement from the Higgs sector. 

We finally studied CP-odd variables at e + e~ linear colliders. We found that the CP decay 
asymmetry constructed from the final state kinematics is typically about 0.1 — 1.5%. It will be 
challenging to experimentally observe this rather small asymmetry. 



Acknowledgments 

We thank the organizers of the "30 Years of Supersymmetry" Workshop in Minneapolis, where 
some of this work was completed. We would like thank Maxim Pospelov for many discussions, 
along with Norval Fortson, Mike Romalis, Jan Kalinowski and Hitoshi Murayama. TF would also 
like to thank the TPI in Minneapolis for its gracious hospitality this fall. This work was supported 
in part by DOE grant DE-FG02-95ER-40896, in part by the University of Wisconsin Research 
Committee with funds granted by the Wisconsin Alumni Research Foundation, and in part by 
DOE grant DOE-EY-76-02-3071. 

References 

[1] T. Falk and K. Olive, Phys. Lett. B375, 196 (1996). 



26 



[2] T. Ibrahim and P. Nath, Phys. Lett. B418, 98 (1998); T. Ibrahim and P. Nath, Phys. Rev. 
D57, 478 (1998); Erratum-ibid. D58, 019901 (1998); Erratum-ibid. D60, 079903 (1999); 
Erratum-ibid. D60, 119901 (1999). 

[3] T. Ibrahim and P. Nath, Phys. Rev. D58, 111301 (1998); Erratum-ibid. D60, 099902 (1999); 
M. Brhlik, G. J. Good and G.L. Kane, Phys. Rev. D59, 115004 (1999). 

[4] V. Barger, T. Han, T. Li and T. Plehn, Phys. Lett. B475, 342 (2000). 

[5] S. Y. Choi, A. Djouadi, H. S. Song and P. M. Zerwas, Eur. Phys. J. C8, 669 (1999); J. L. Kneur 
and G. Moultaka, Phys. Rev. D61, 095003 (2000); S. Y. Choi, A. Djouadi, M. Guchait, 
J. Kalinowski, H. S. Song and P. M. Zerwas, Eur. Phys. J. C14, 535 (2000). 



[6] S. Abel et al. [SUGRA Working Group Collaboration], Eiep-ph/0003154 



[7] J. Ellis, J. S. Hagelin, D. V. Nanopoulos and M. Srednicki, Phys. Lett. B127, 233 (1983); 
V. Barger, R. W. Robinett, W. Y. Keung and R. J. Phillips, Phys. Lett. B131, 372 (1983); 
S. Dawson, E. Eichten and C. Quigg, Phys. Rev. D31, 1581 (1985); W. Beenakker, M. Klasen, 
M. Kramer, T. Plehn, M. Spira and P. M. Zerwas, Phys. Rev. Lett. 83, 3780 (1999). 

[8] T. Tsukamoto, K. Fujii, H. Murayama, M. Yamaguchi and Y. Okada, Phys. Rev. D51, 3153 
(1995); J. L. Feng, M. E. Peskin, H. Murayama and X. Tata, Phys. Rev. D52, 1418 (1995); 
S. Y. Choi, A. Djouadi, H. Dreiner, J. Kalinowski and P. M. Zerwas, Eur. Phys. J. C7, 
123 (1999); G. Moortgat-Pick, H. Fraas, A. Bartl and W. Majerotto, Eur. Phys. J. C9, 521 
(1999). 

[9] M. V. Romalis, W. C. Griffith, and E. N. Fortson, |iep-ex/001200T| ; see also J.P Jacobs et 
al, Phys. Rev. Lett. 71, 3782 (1993) and J.P. Jacobs et al, Phys. Rev. A52, 3521 (1995). 

[10] S. Dimopoulos and D. Sutter, Nucl. Phys. B452 (1995) 496. 

[11] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 

[12] I. Bigi and N.G. Uraltsev, Sov. Phys. JETP 100 (1991) 198; M. Pospelov, Phys. Rev. D58, 
097703 (1998). 

[13] J. Ellis, S. Ferrara, and D.V. Nanopoulos, Phys. Lett. 114B, 231 (1982); W. Buchmiiller and 
D. Wyler, Phys. Lett. 121B, 321 (1983); J. Polchinski and M. Wise, Phys. Lett. 125B, 393 
(1983); F. del Aguila, M. Gavela, J. Grifols, and A. Mendez, Phys. Lett. 126B, 71 (1983); 
D.V. Nanopoulos and M. Srednicki, Phys. Lett. 128B, 61 (1983); M. Dugan, B. Grinstein and 



27 



L. Hall, Nucl. Phys. B255, 413 (1985); C.-S. Huang and W. Liao, Phys. Rev. D62, 016008 
(2000). 

[14] Y. Kizukuri & N. Oshimo, Phys. Rev. D45, 1806 (1992); D46, 3025 (1992). 

[15] T. Falk, K. A. Olive, M. Pospelov and R. Roiban, Nucl. Phys. B560, 3 (1999). 

[16] K.F. Smith et ai, Phys. Lett. B234, 191 (1990); I.S. Altarev et ai, Phys. Lett. B276, 242 
(1992). 

[17] E.D. Commins et ai, Phys. Rev. A50, 2960 (1994). 

[18] T. Falk and K.A. Olive, Phys. Lett. B439, 71 (1998). 

[19] T. Falk, K.A. Olive and M. Srednicki, Phys.Lett. B354, 99 (1995). 

[20] M. Pospelov and A. Ritz, Nucl. Phys. B573, 177 (2000); M. Pospelov and A. Ritz, Phys. 
Rev. Lett. 83, 2526 (1999). 

[21] M. Pospelov and A. Ritz, |hep-ph/0010037 . 

[22] S. Weinberg, Phys. Rev. Lett. 63, 2333 (1989). 

[23] A. Pilaftsis, Phys.Lett. B471, 174 (1999). 

[24] T. Ibrahim and P. Nath, Phys. Rev. D61, 093004 (2000); M. Brhlik, L. Everett, G. L. Kane 
and J. Lykken, Phys. Rev. Lett. 83, 2124 (1999). 

[25] K. Inoue, A. Kakuto, H. Komatsu and S. Takeshita, Prog. Th. Phys. 68, 927 (1982); M. Drees 
and M. Nojiri, Phys. Rev. D47 376, (1993); W. de Boer, R. Ehret, W. Oberschulte, and D.I. 
Kazakov, Z. Phys. C67, 647 (1995); S.P Martin and M.T. Vaughn, Phys. Rev. D50, 2282 
(1994). 

[26] R. Garisto and J. D. Wells, Phys. Rev. D55, 1611 (1997). 

[27] E. Accomando, R. Arnowitt and B. Dutta, Phys. Rev. D61, 115003 (2000). 

[28] Recent LEP limits on SUSY particles and Higgs bosons can be found in the November LEPC 
talks, which can be found at 

http : //delphiwww. cern. ch/~of f line/physics_links/lepc .html 



[29] S. Heinemeyer, W. Hollik and G. Weiglein, Phys. Rev. D58, 091701 (1998); see also 
R.-J. Zhang, Phys.Lett. B447, 89 (1999). 



28 



[30] S.Y. Choi, H.S. Song and W.Y. Song, Phys. Rev. D61, 075004 (2000). 

[31] V.M. Khatsimovsky, LB. Khriplovich and A.S. Yelkhovsky, Ann. Phys. 186, 1 (1988). 

[32] O.P Sushkov, V.V. Flambaum and LB. Khriplovich, Zh. Eksp. Theor. Fiz. 87, 1521 (1984) 
(Sov. Phys. JETP 60, 873 (1984)). 

[33] D. A. Dicus, S. Nandi, W. W. Repko and X. Tata, Phys. Rev. Lett. 51, 1030 (1983); A. Bartl, 
H. Fraas and W. Majerotto, Z. Phys. C30, 441 (1986); A. Bartl, H. Fraas and W. Majerotto, 
Nucl. Phys. B278, 1 (1986); E. Accomando et al. Phys. Rept. 299, 1 (1998); T. Li, Ph.D. 
thesis (2000). 

[34] H. Martyn and G. A. Blair, 4th International Workshop on Linear Colliders, Sitges, 1999. 
[hep-ph/9910416 . 



[35] G. Moortgat-Pick and H. Fraas, Acta Phys. Polon. B30, 1999 (1999). 

[36] J.L. Feng and T. Moroi, Phys. Rev. D56, 5962 (1997); V. Barger, T. Han and J. Jiang, 
|hep-ph/0006223 . 



[37] J. Kalinowski, private communication. 

[38] Eugene Commins, private communication. 

D. DeMille et al, "Investigation of Electronic Transitions in PbO", submitted to Phys. 
Rev. A. 



[40] D. A. Demir and M. B. Voloshin, |iep-ph/0012"T23 



29 



